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Abstract 

We show how to extrapolate an analytic function (or a smooth signal) by multiplying and dividing 
its values on geometric sequences that collapse to a point. 
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1. Introduction 

Digital signal processing normally uses signal samples that are uniformly spaced in time; the 
sample measurement times form an arithmetic sequence [3]. This article shows how to extrap- 
olate an analytic function (or a smooth signal), by multiplying and dividing function values (or 
signal samples) at points that form sets of geometric sequences. 

The starting point for this work was a trigonometric identity attributed to Euler, written in a 
particular way: 

00 1 r* x 

n cos (l) = *j cos(s)ds - 

Eqation 11.11 possesses a certain parallelism. In forming the simple integral on the right, 
the cosine function is sampled on a dense arithmetic sequence, the values are added, and the 
sum is normalized. On the left, the cosine function (already normalized to unity at the origin) 
is sampled on a geometric sequence, and the values are multiplied. The equation also has a 
predictive capability: on the right we have the average value of cosine on the interval [0, x], 
whereas on the left we evaluate cosine only in the first half of the interval. Lastly, since the 
equation has only one function present, one might wonder whether there is a generalization for 
other functions. This prompted us to look for a similar equation, which we hoped would predict 
the value of a function at one point by multiplying prior geometric values. 

The solution we found does indeed multiply function values from geometric sampling. But 
the function is evaluated on multiple geometric sequences. And, strengthening the parallelism, 
the values are multiplied as the geometric sequences become dense by collapsing to a point. 
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2. An Extrapolation Tool 



We call our result the Geometric Multiproduct. 

We let x represent the point to which we want to extrapolate the function /. To assure 
convergence of infinite products [4], we require that / be normalized to unity at the origin and 
that the geometric ratio (1/r) be less than one. Our result is expressed in the following theorem. 

Theorem 2.1. Suppose that f: R — > R is analytic and nonzero on the interval [0, x], that /(0) 
I and that r > 1. Let S be a set of positive integers, S c N*. If we define 



x„ (S, r, x) 



i/k 



keS 



X 



and P f (S,r,x) = f\ [f(x n (S,r,x))] 



(\S\-\) 



n=\S\ 



then f{x) = lim 



II Pf(.S,r,x) 

\S\ odd 



^ n P f (S,r,x) 

\S\ even 



(2.1) 



(2.2) 



(2.3) 



Discussion 



Equation l2.1l defines the geometric sequences {x n }°°_ v There is a geometric sequence {x n } 
for each finite subset S of the natural numbers N* (positive integers). The points {x„} 
depend upon x and the geometric ratio r > 1, but also have a coefficient that depends upon 
the integers in the set S. For simplicity, we have assumed that the geometric sequences 
converge to the origin. 



Equation |2j2] uses the function / to calculate a product Pf for each set S. Pf is obtained 
by multiplying values of the function /, evaluated on the geometric sequence for the set S. 
The index n starts at \S\, not 1, and each factor of / has a multiplicity These details 

will be explained later. 



• Equation 12 . 3 1 calculates the value of the function / at the point x by multiplying and di- 
viding factors of Pf, There is one factor Pf in the numerator for each subset S of N* 
containing an odd number of integers, and one factor Pf in the denominator for each sub- 
set S of N* containing an even number of integers. After multiplying and dividing factors 
of Pf, we take the limit as r approaches 1 from above, r I 1. 

Notice that all points {x n } in all geometric sequences aproach zero as r I 1. We shall prove 
that, although /(0) = 1, the formula produces f(x), not 1. 



3. Two Simple Examples 

Before proving our result, we show how to use the Geometric Multiproduct by extrapolating 
two simple functions. To make the calculations manageable, we must truncate the infnite prod- 
ucts and pick a value for r near 1 . A truncated Geometric Multiproduct gives us an estimate for 
the value of a function. 
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Figure 1: Geometric Multiroducpt for cos(:t). Solid line is the estimate; dashed line is the function. The estimate has 
Smax = (2,4), n max = 10, and r = VI 
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Figure 2: Geometric Multiproduct estimate for 1 + ^sin(x). Solid line is the estimate; dashed line is the function. The 
estimate \msS max = (1,2,3,4), n max = 40, and r = 2. 



We first perform an extrapolation for cos(x). For this example, we let the index n run only 
from 1 to n max = 10. We truncate the sets \S) at S max = {2,4}. This means we use the power set 
of S max '■ {2}, {4}, and {2, 4). (Later we will show why we only need to use even integers for even 
functions.) This gives 



rC =1 cos [(r 4 -l)^x/r"] 



cos(x) * hm !P cos ({2}, r, x) ■ - — — — 

f , cos ({2,4},r,x) 

(- 10 

lim iflcoskr 2 - l) 1/2 x/r"l 

~l + \U 1 J J Y\nLi (cos [(r 2 - ly/Hr 4 - D 1 ' 4 x/f]} ' 

For accuracy, we should truncate the geometric sequences of cosine arguments near their 
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limit point of zero. This is satisfied if r w is large. To get r w — 2 5 , we set r = V2. Graphing both 
sides, we obtain the result shown in Fig. Q] The agreement between the Geometric Multiproduct 
estimate and cos(x) is striking, given how drastically we truncated the infinite products and how 
far r — y/2 is from 1. 

Figure [2] shows an extrapolation for the function 1 + 5 sin(x). If we wanted to extrapolate 
sin(x), we might use this function instead, to assure that /(0) = 1 and / > 0. The extrapolation 
has S max = {1,2,3,4}, n max = 40, and r = 2. 

Increasing S max would give more oscillations to the estimate for this function. In this esti- 
mate, calculation of / at one point involves multiplication of more than 75,000 factors of earlier 
values of /. We see again that Geometric Multiproducts are highly tolerant of r-values that are 
quite far from 1, at least in percentage terms. 



4. Proof of the Geometric Multiproduct 

Proof. We divide the proof of Theorem l2.1l into five short lemmas. 



4.1. An Invariance 

We first define a geometric sequence and a corresponding function that have a very useful 
property. 

Lemma 4.1. Let {x n }™ =x be the geometric sequence x n (k,r) — (r - l) 1 ^ x/r n , with r > 1. Let 
fic — exp^jc*), where q is a constant and k is an integer. Then 



\~[fk(x„) = fk(x). 



(4.1) 



Here 1 /r is the geometric ratio, and x is the point at which we want to calculate the value of 
/. Loosely speaking, the function fk(x) is invariant under multiplication of function values on 
the sequence x„(k, r). 



Proof. The proof is straightforward algebra: 



Ik* — ? — =1 l exp 

n=i x ' n=\ 



Ck 



(r* - l) 1 ^ 



exp 



c k x 







)" 


= exp 






f 


= exp 



11=1 * 
QXV-D- 11,1 



1 - 1/r* 



exp(Q^) = f k (x). 



□ 



Since fk{x) is independent of r, we may take the limit r I 1 on the left-hand side of the 
equation. 

" l( r k _ \)Vk x \ 



n=\ 

The limit r I 1 makes the geometric sequence dense, because as r I 1 all points of the geometric 
sequence approach the origin: lim(r k - l) l ^ k x/r n = 0. Eq. l4.2l savs we can multiply values of /j 

evaluated on a geometric sequence within an infinitesimal neighborhood of zero and, although 
fii(0) = 1, we do not get 1. We get the value of fa at the point x. 
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4.2. Component Factors of a Function 

We write a normalized function / as a product of /i's. If /: C — > C is holomorphic and 
nonzero in some simply-connected open neighborhood, therein log(/) is also holomorphic and 
can be written as a power series (!]]: 



f(x) = e hg(f > = exp(^ c k x k ) = \\ f k (x), 



(4.3) 



k=\ 



where, as before, fk(x) = exp(Qx*) with k an integer. By restriction, the statement is also true 
for a positive real-analytic function / on open interval of the real axis. 

We call the factor fk(x) the k-th component of the function/. For an even function, only even 
values of the index k would appear. (This is why only even integers occur in the sets {S} for an 
even function.) 

It would be convenient if we could produce fa merely by geometrically sampling / on the 
order-£ sequence. However, when we sample /, we sample not only its fa component, but also 
all of the other components of /. So we need to understand what happens if we sample the j-th 
component of / on the order-fe sequence. 

Lemma 4.2. Let fj(x) = exp(c ; x j ). Then 



8? I I 1 



(r* - 1) 



l Ik 



fk(x) ifj = k 

1 if j > k or cj — 

if j < k and Cj < 

oo if j < k and cj > 



(4.4) 



Proof. 



n*(^)-n4( 

H-1 \ f .1-1 \ 



>(r k - l) l/k x\ j 



= exp 



= exp 



n=i * ' 



CjX 



j(r k -iy /k 
ri- 1 



Letting r — 1 + 6, using the binomial theorem, keeping lowest order terms in 6, and letting r I 1 
(i.e., S i 0), we obtain 



(r* - iy /k . 



[(l+flt-l]j/* 

hm exp \CjX J : 



I , [1 +k6 + 0(6 2 )- \V lk 

= hm exp c , x J r 

sio F \ 1 (1 + j5 + 0(6 2 ) - 1) 

= lim exp \ciX j ■ ■ lim S (j/kyl \ 

sio \ j 610 ) 



The cases in Eq. l4.4l follow directly. 



□ 



4.3. A Key Intermediate Result 

According to Eq. 14.41 taking r I 1 (i.e., 6 I 0) eliminates the effect of components fj with 
j > k from the order-£ geometric sampling of /: they each provide a factor of 1 . The effect of 
the finitely-many components with j < k can be divided out before the limit is taken. 

Lemma 4.3. Let fa(x) = exp(cjX k ). Then 



/*(*) 



«=i 



(r* - 1) 



(r* - 1) 



l/k 



(4.5) 



Proof. To express /j in terms of /, we first augment fa with the other components necessary to 
construct /. 



/*(*) 



CO . 

M-1 ,>(■ V 



(R k - iy /k x 

~R" 

(R k - iy /k x 



ft 



(R k - iy /k x 



ok 



(r* - l) 1/k 



in* 

i<k 



(R k - iy /k 

R" 
(4.6) 



This equation is true for any value of R. The last factor in brackets in the numerator is unity. If 
we set R = r, we can take the limit of the entire righthand side. This gives Eq. I4.5I □ 



Note that, as empasized in equation I4.61 it is only the higher-order components (j > k) 
that require us to take the limit r I 1. Since we expect the high-order components to be less 
significant than low-order ones, we might also expect our results to depend only weakly upon 
taking the limit r { I. Numerical examples show that we do not need to use an r-value very close 
to 1. 



4.4. Calculating the Components of a Function 

We want to express the value of / at x solely in terms of earlier values of / itself. We can 
calculate f(x) if we can calculate its components fk(x). To calculate the denominator of fa in 
Eq. ( 14.5b . we must know how to calculate f for all i < k. But where does f come from? Answer: 
the formula for fa can be used recursively to calculate its own denominator. 

The first component has no denominator and can be calculated solely from /: 

oo 

/!(*) = lim U f 

Hi =1 

The second component fa will have / in its numerator and fa in its denominator, but fa can be 
expressed in terms of /: 



oo 

f 2 (x) = urn n / 

n\ 1 l . 

« 2 = 1 



(r- 1) Vr 2 - 1 x 
6 



(r-l)x 



Vr 2 - 1 x 



fi'i 



lim 

rj.1 



Ylf 

m .n-> — 1 



Vr 2 - 1 x 



If 



The number of factors rapidly increases: f$ will have fx and f\ in its denominator, but the f in 
the denominator will spawn an f\ in the denominator of the denominator (i.e., in the numerator). 
After expressing each f\ in terms of /, 



Mx) - lim (1 / 



r" 3 



// 



(r- l)(r 2 - l) 1/2 (r 3 - 1) 



1/3 



(r- l)(r 3 - 1) i/3 jc 



(r 2 



1} l/2 (r 3 



1) 



1/3 



r « 2 4 



There is a shortcut notation that can help us cope with the complexity. The third component 
can be expressed as 

,. n /({3}°x)/({l,2,3}ox) 

fifjt) = lim , 

Mi 1 1 /({l,3}ox)/({2,3}ox)' 

in which {. . . , ./',...} ° x indicates the presence of a factor (V - l) 1 ^/'""' multiplying x, and in 
the product Y\ all indices nj run from 1 to oo. 
As k gets larger, we begin to see a pattern. 

• In the numerator of fk we have all sets of positive integers whose greatest element is k and 
have an odd number of elements. 

• In the denominator we have all sets of positive integers whose greatest element is k and 
have an even number of elements. 

We prove that this pattern always holds. 

Lemma 4.4. Let the index set [St) be the family of sets of positive integers, each of whose largest 
element is k: Su C N* and max^t) = k. Then, with the previous assumptions and notation, 



CO oo 

Mx) =lim f| ]~[ f(S k o X )/ p| p[ f(S k ox) 



(4.7) 



[5*1 odd 



flj, . . 



15* I even 



■ , flj, . . 



Proof. The proof uses induction on k. We have shown that the assertion holds for small k. We 
show that if it holds for all i < k, it also holds for k. Eq. l4.5l can be written 

f({k}o X ) 



" k - 1 i<k 



({k} o x) 



n n n 

i<k \Sj\ even «^=1 



= ]imf]f((k}°x) ■ 

ri 1 1 - 1 - 

n k =\ 



n 

Hj, . . . 



f([SiU{k)]ox) 



n n n 

i<k |5,| odd «t=l 



n 



/(KuR)]oi) 



where we have used Eq. l4.7l to express f, for all i < k. In the denominator, after unioning each 
odd cardinality set Si , i < k with the set {k}, we have all of the even cardinality sets in {St} in the 
denominator. Similarly, we end up with all of the odd cardinality sets in {Sk} in the numerator, 
including the singleton set {k}, which appears in the first factor. This gives Eq. 14.71 completing 
the induction. □ 
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4.5. Simplification: Consolidating Exponents 

We can simplify Eq. l4.7l bv consolidating all indices nj, j e Sk, into a single index n. 

Lemma 4.5. Let Sk c N* with max(Sk) — k, where k is a positive integer. With the previous 
assumptions and notation, each component fk(x) off is given by 



Ft Pf(S k ,r,x) 

, .. \Sk\odd 

fk(x) = hm , 

M n r Pf{S k ,r,x) 

\Sk\ even 



(4.8) 



where V /(Sk, r, x) is defined by Eauations \2.2\ and \2. 1 
Proof. If we express Eq. l4~7l explicitlv. 



fk(x) = lim 

rj.1 



n n /(|n,.., 

|S S | odd ..,n/,...= l 

~n n 

|*Sjfc| even rij,,..= 1 



(4.9) 



In both the numerator and denominator of fk(x), there is a product 



n / 



...,»„...= i 



- n 



B=]Stl 



n 



n (^ - 1) 



n / 



...,»,,... = ! 



fl (rJ_1) 



i/; 



Z »j 
r j e S t 



Uj 



= f[ 7M..SV'-..v) . 



n=l-Si 



The binomial coefficient (j^T^i) is the number of ways in which the \Sk\ indices {. . . , tij, 
Sk}, can add up to n (Ref. [2]). The index n — Yu n j starts at |»Sjtl because each rij is at least 1, 
and there are \Sk\ of them. Substituting this result in both numerator and denominator of Eq. |4.9 
gives Eq.l4~8l □ 



To obtain fix), we multiply all fk(x) for k e W. This effectively removes the restriction 
max(Sk) = k in Eq. 14.81 giving us the Geometric Multiproduct for fix), Eq. 12.31 This ends the 
proof of Theorem l2.ll □ 

Mathematically, the Geometric Multiproduct does what Taylor series do. It uses information 
about a function from an infinitesimally small neighborhood of one point to calculate the value 
of the function at another point. 



4.6. Practical Applications 

In engineering applications, to sample a signal means to measure its magnitude periodically, 
e.g., at moments that constitute an arithmetic sequence in time, or at uniformly separated places. 
The signal is often a voltage that represents a physical quantity, such as the pressure of a sound 
wave or the intensity of light at one point. Periodic sampling of a signal is widely used in many 
modern inventions, such as digital telephone networks and digital cameras. The samples are 
quantized and digitized (rounded to an integer and converted to a binary number) to facilitate 
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computation and transmission. Subsequently, the digitized samples can be converted back to an 
analog signal, e.g., speech or image. In finance, recording the daily closing price of a stock is 
another example of sampling. 

Sampling has an important advantage over analytic methods. If we can do all of our calcula- 
tions with the digitized samples alone, we never need to determine the analytic form of the signal. 
Ref. ||3t] describes how to process a signal by performing calculations with digitized arithmetic 
samples, for example, to remove unwanted frequency components in a digitized speech signal. 

The Geometric Multiproduct shares the same advantage. Even if we don't know the analytic 
form of a signal, we can still make predictions of its future value, by using geometric samples: 
measurements of the signal at points in time that form geometric sequences. For this reason, it 
could potentially have many applications in engineering, science, finance and any other field in 
which forecasting plays an important role. 
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